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We present a neutron scattering investigation of Ce\- X Y x h\s as a function of chemical pressure, 
which induces a transition from heavy fermion behavior in CeAl3 (Tk = 5 K) to a mixed valence state 
at x = 0.5 (Tk = 150 K). The crossover can be modeled accurately on an absolute intensity scale by 
an increase in the k-f hybridization, Vkf, within the Anderson Impurity Model. Surprisingly, the 
principal effect of the increasing Vfe / is not to broaden the low-energy components of the dynamic 
magnetic susceptibility but to transfer spectral weight to high energy. 

PACS numbers: 71.27.+a,71.70.Ch,75.40.Gb,78.70.Nx 



The Anderson Impurity Model (AIM) has been invoked 
to describe the thermodynamic, transport, and spectro- 
scopic properties of strongly correlated electron systems 
for nearly fifty years [1]. It is believed to contain the 
essential physics of materials in which a narrow band of 
quasi-localized electrons are hybridized with a broader 
band of itinerant electrons, and provides a mechanism 
describing the formation and screening of local moments 
in metallic systems. The AIM has been most widely ap- 
plied to describe fluctuating moments in rare earth f- 
electron systems, ranging from strongly hybridized mixed 
valent systems to more weakly hybridized heavy fermion 
systems [2]. Although it is a single- impurity model, it 
has been successfully applied to concentrated rare earth 
systems, except at the lowest temperatures where lattice 
coherence cannot be ignored. 

In spite of its widespread use in heavy fermion and 
mixed valence physics, there have been no direct quan- 
titative comparisons of ab initio calculations of the dy- 
namic magnetic susceptibility, Imx(w), with experiment 
covering the crossover between the two regimes. In re- 
cent years, it has become possible to perform detailed 
calculations of Imx(w) using the Non-Crossing Approx- 
imation (NCA) in the presence of crystal field splittings 
of the 4/-electron ground multiplet [3, 4], but these theo- 
retical advances have never been tested against inelastic 
neutron scattering data, which directly measure Imx(w). 
The purpose of this letter is to present the first detailed 
quantitative comparison of AIM/NCA calculations with 
neutron scattering data as a function of the hybridization 
strength, Vkf, covering the transition from heavy fermion 
to mixed valence behavior. There is excellent agreement 
on an absolute intensity scale between theory and exper- 
iment, although the relative strength of the transverse 
and longitudinal susceptibility requires an adjustment for 



anisotropic hybridization. The low-frequency broadening 
of the magnetic response is typically proportional to the 
Kondo temperature, Tk, so it is surprising that increas- 
ing the hybridization does not significantly increase the 
width of the low-energy contributions to Imx(w). In- 
stead, the transition is predominantly characterized by a 
transfer of spectral weight from the two low-energy com- 
ponents, respectively a narrow quasielastic response and 
a crystal field excitation, to a broad high-energy tail. 

The variation in hybridization is achieved through the 
application of chemical pressure on the canonical heavy 
fermion system CeAl 3 [5]. Doping lanthanum onto the 
cerium sites expands the lattice, progressively localizing 
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FIG. 1: The unit cell volume, the static susceptibility mea- 
sured at 8K and the magnetization in a field of 12T measured 
at 2K of Cei^yLa^Als and Cei-^Y^Ala. The solid circles 
are measurements and the open squares are predictions of a 
localized crystal field model. The lines are guides to the eye. 
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FIG. 2: S(Q,w) of Ce\~ x Y x A\3 as a function of x measured at an average scattering angle of (a) 60° and (b) 15°, with an 
incident energy of (a) 3.12 meV and (b) 35 meV. The solid lines are fits to a quasielastic Lorentzian lineshape (dashed line) 
and an elastic peak (dotted line) convolved with the instrumental resolution (FWHM ~ 90/ieV). The inset shows the fitted 
values of the static susceptibility, xo, and the half-width of the quasielastic peak, V. 



the 4/ electrons while doping with yttrium compresses 
the lattice, inducing a transition from the heavy fermion 
behavior of pure CeAl 3 , characterized by a Kondo tem- 
perature of Tk = 5K, to a mixed valence state in 
Ceo.5Yo.5Al3, where Tk = 150K. Comparison with the 
cell volumes of pure CeAl 3 under pressure shows that 
the chemical pressure in Cei_ x Y x Ai3 at x = 0.5 is equiv- 
alent to 27 kbar. The bulk measurements in Figure 1 
illustrate how chemical pressure, and the consequent in- 
crease in hybridization, reduces the static susceptibility 
and high-field magnetization from the values predicted 
in a pure crystal field model [6]. 

Polycrystalline samples of Cei-^Y^ALj, with x = 0.0, 
0.05, 0.1, 0.2, 0.3, 0.4, and 0.5, Ce^La^ with y = 0.2, 
0.3, 0.4, 0.5, and 0.95, along with non-magnetic samples 
of LaAl3 and Lao.5Yo.5Al3, were made by arc-melting sto- 
ichiometric quantities of the constituent elements, with 
a total mass of approximately 30 g for each sample, fol- 
lowed by annealing in a vacuum at 900° C for four to five 
weeks. Powder neutron diffraction confirmed that each 
sample was single-phase. The inelastic neutron scatter- 
ing experiments were performed on the IN6 spectrometer 
at the Institut Laue Langevin, using an incident neutron 
energy of 3.12 meV, and the LRMECS spectrometer at 
the Intense Pulsed Neutron Source, in Argonne National 
Laboratory, using incident energies of 35 and 80 meV. All 
spectra were normalized on an absolute intensity scale 
using a standard vanadium plate and corrected for self- 
shielding, absorption, and form factor. Our own compar- 
ison with inelastic scattering standards has shown that 
it is possible to achieve accuracies of 1% or better, a pre- 
cision that is confirmed here in our comparisons of data 



taken at different incident neutron energies. 

Figure 2 shows corrected data from IN6 and LRMECS 
at a temperature of 8 K covering the low and intermedi- 
ate energy range. In Figure 2(a), the magnetic response 
is fitted to a single quasielastic Lorentzian lineshape con- 
volved with the instrumental resolution. 

S(Q,u) cx F 2 (Q) [„( W ) + 1] gp^y (1) 

where F(Q) is the Cc 3+ magnetic form factor, n(u>) is the 
Bose population factor, \o is the static susceptibility, and 
T the half width. 

With increasing x, there is only a modest increase in 
the linewidth V from 1.06 meV at x = to 1.26 meV at 
x = 0.3, but a dramatic reduction in the spectral weight, 
XOj becoming negligibly small at x = 0.5. The energy 
range of Fig. 2(b) covers the crystal field excitation, 
which is observed at 6.7 meV in CeAl 3 . With increas- 
ing hybridization, the crystal field peak falls sharply in 
intensity, with a slight shift to higher energy, and is dif- 
ficult to resolve at x > 0.3. 

At all values of x, there is an additional component 
to the magnetic response in the form of a high-energy 
tail, extending to greater than 70 meV, that cannot be 
accounted for by the low-energy quasielastic Lorentzian 
or the crystal field excitation. It is also just visible at 
y = 0.3, but it has negligible intensity at y = 0.4. At 
these energies, multiple scattering produces a significant 
phonon contribution to the measured inelastic scattering 
and must be subtracted before comparing to theory. Fig- 
ure 3 illustrates this with data at x — 0.4, which shows 
significant additional magnetic intensity when compared 
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to the non-magnetic LaAl3 data. The inset to Figure 3 
shows that the high-energy magnetic component, inte- 
grated from 40 to 70 meV, falls with increasing cell vol- 
ume, i.e., with decreasing hybridization, becoming neg- 
ligible in Ceo.6Lao.4Al3. 

Qualitatively, the observations are consistent with 
AIM/NCA calculations, which follow the method of Ku- 
ramoto [4] and Bickers et al. [3] (Figure 4). We as- 
sume a Gaussian conduction band of half-width W = 3 
eV, as in Ref. [7] with the ground state of the Ce 3+ 
4/ electrons at Ef = —2 eV. The 14- fold degeneracy of 
the /-states is lifted by the free-atom spin-orbit coupling 
and the crystal field potential, which, in hexagonal point 
group symmetry, is diagonal in J z . The J = | multiplet 
is therefore split into three Kramers doublets, IV = (±5), 
Tg = |±|), and Tg = |±|). In an earlier publication, 
we showed that the Tg doublet is the ground state, with 
the other two doublets nearly degenerate at an energy of 
approximately 7 meV [6]; in the calculation, the T$ and 
T 7 doublets were assumed to be at 6 meV and 7.5 meV, 
respectively. The hybridization matrix element was as- 
sumed to be isotropic and its value was treated as a free 
parameter in order to simulate the effect of increasing x. 

For low values of the hybridization, the AIM/NCA cal- 
culations show a sharp quasielastic response, well-defined 
crystal field excitation, and a broad high energy tail. 
The high-energy tail in the spectrum is a consequence 
of the dynamic screening of the local moments by the 
conduction electrons [8] . The screening of the magnetic 
moments can be considered as an iterated orthogonality 
catastrophe [9], which occurs since the conduction elec- 
trons experience a sequence of local spin-dependent po- 
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FIG. 3: S(Q,o;) of Ceo.6Yo.4Al3 (open circles), Ceo.6Lao.4Al3 
(solid circles) , and LaAl3 (diamonds) measured at an average 
scattering angle of 15° with an incident energy of 35 meV 
(left axis) and 80 meV (right axis) at 8 K, normalized on an 
absolute scale. The shaded area shows the estimated mag- 
netic scattering. The inset shows the magnetic scattering in- 
tegrated from 40 to 70 meV as a function of La and Y doping. 
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FIG. 4: Theoretical calculations of S(oj) determined by the 
Anderson Impurity Model using the Non-Crossing Approxi- 
mation at a temperature of 8 K. 

tcntials caused by successive flips of the local spin. This 
process results in a power-law decay of the spectrum at 
high energies, similar to that found in the x-ray edge 
problem [10]. 

As the hybridization increases, the calculated 
quasielastic and crystal field contributions gradually de- 
crease in intensity, but their spectral weight is transferred 
into the high energy tail. This agrees with the experimen- 
tal trends with increasing x. Since the theoretical predic- 
tions cannot be described by simple analytic functions, 
we have put this comparison on a quantitative footing 
by directly comparing integrals over energy ranges corre- 
sponding to the three components in the scattering with 
both the experimental and calculated intensities placed 
on an absolute scale with no adjustable parameters. The 
quasielastic component is represented by an integral from 
-3 to 2 meV (using the fitted profile to exclude the elastic 
nuclear scattering) , the crystal field component from 5 to 
15 meV, and the high energy component from 40 to 70 
meV. For the latter two components, similar integrals 
were performed on LaAl 3 and Lao.5Yo.5Al3, and sub- 
tracted after interpolating for the correct yttrium con- 
centration. The integrals in the high-energy tail are not 
affected by this interpolation since the phonon scattering 
at these energies is dominated by the aluminum contri- 
bution. 

Figure 5(a) shows that two contributions to Im%(w), 
the narrow quasielastic and high energy components, are 
in excellent agreement with the AIM/NCA calculations 
with V kf = 443 meV at x — and 477 meV at x = 0.5. 
It appears that there is a nearly linear relation between 
Vkf and x, although there is no a priori reason for this 
to be so. The important result is that both the low en- 
ergy and high energy regions of scattering have absolute 
cross sections that are accurately predicted at each value 
of Vkf over the entire transition from heavy fermion to 
mixed valence states. 
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FIG. 5: The magnetic cross section of Cei-^Y^Ala as a func- 
tion of x integrated over three energy regions, -3 to 2 meV 
(triangles), 5 to 15 meV (circles) and 40 to 70 meV (squares), 
and corrected for the form factor. The solid lines are the 
same integrals determined from the theoretical calculations 
using the AIM/NCA model as a function of Vkf, using a lin- 
ear mapping of Vt,/ to x. The dashed line is the same as the 
solid line shifted to illustrate the effect of increasing Vfc/ by 
~ 20 meV at each value of x. 



ity. There is some evidence that anisotropic hybridization 
may also be important in other strongly correlated elec- 
tron systems, such as CeCug-zAu^lS], where the mag- 
netic fluctuations appear two-dimensional, and URu 2 Si2 
[14], where we have proposed that the AKM plays an 
important role in generating hidden order [12]. 

In conclusion, a quantitative comparison of inelastic 
neutron scattering data with NC A calculations of the An- 
derson Impurity Model have provided, for the first time, 
a consistent description of the evolution of the spin dy- 
namics from heavy fermion to mixed valence behavior. 
Because of the non-analytic character of the Kondo cou- 
pling, a 10% increase in the k-f hybridization produces an 
increase in the energy scale of spin fluctuations of more 
than an order-of-magnitude. However, the most signifi- 
cant result is that the transition between the two regimes 
does not primarily result from a general broadening of the 
dynamic magnetic susceptibility but by a significant shift 
of spectral weight from from the narrow components that 
characterize well-localized /-electron states to a broad 
magnetic response that is present at all hybridization 
strengths, including heavy fermion compounds such as 
CeAl 3 . 

This research was supported by the U.S. Department 
of Energy, Office of Basic Energy Sciences, Division of 
Materials Sciences and Engineering under Award # DE- 
AC02-06CH11357. 



On the other hand, Figure 5(b) shows that the 
measured intensities in the intermediate energy region, 
which, at lower values of x, includes the crystal field ex- 
citation, are much lower than predicted if we use the 
same values of Vkf- Reasonable agreement between the- 
ory and experiment can only be achieved in this energy 
range if we increase the value of / by approximately 
20 meV at each value of x. To explain this apparent 
discrepancy, it is important to know that, because the 
crystal field potential is diagonal in J z , the quasielastic 
response is purely longitudinal, whereas the crystal field 
excitation is purely transverse in character. The anomaly 
in the crystal field intensities therefore implies that the 
k-f hybridization in CeAl3 is anisotropic. 

In earlier investigations, we used the Anisotropic 
Kondo Model (AKM) to explain the anomalous spin dy- 
namics in Cei-^La^Als [11, 12]. However, the present re- 
sults provide the first direct evidence of the anisotropy in 
Vfc / and its impact on the dynamic magnetic susceptibil- 



* Electronic address: R0sborn@anl.gov 
[1] P. W. Anderson, Phys Rev 124, 41 (1961). 
[2] P. Fulde, J Phys F 18, 601 (1988). 

[3] N. E. Bickers, D. Cox, and J. W. Wilkins, Phys. Rev. B 

36, 2036 (1987). 
[4] Y. Kuramoto, Z Phys B 53, 37 (1983). 
[5] K. Andres, J. E. Graebner, and H. R. Ott, Phys Rev Lett 

35, 1779 (1975). 
[6] E. Goremychkin, R. Osborn, and I. Sashin, J Appl Phys 

85, 6046 (1999). 
[7] P. Riseborough, Phys Rev B 67, 045102 (2003). 
[8] P. W. Anderson, Phys Rev 164, 352 (1967). 
[9] P. W. Anderson, Phys Rev Lett 18, 1049 (1967). 
[10] G. D. Mahan, Phys Rev 163, 612 (1967). 
[11] E. Goremychkin, R. Osborn, B. Rainford, and A. Murani, 

Phys Rev Lett 84, 2211 (2000). 
[12] E. Goremychkin, R. Osborn, B. Rainford, T. Costi, 

A. Murani, C. Scott, and P. King, Phys Rev Lett 89, 

147201 (2002). 

[13] H. Lohneysen, A. Rosch, M. Vojta, and P. Woffle, Rev 

Mod Phys 79, 1015 (2007). 
[14] S. Elgazzar, J. Rusz, M. Amft, P. M. Oppeneer, and J. A. 

Mydosh, Nat Mater 8, 337 (2009). 



